Non-linear quantum critical transport and the Schwinger Mechanism 
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Scaling arguments imply that quantum critical points exhibit universal non-linear responses to 
external probes. We investigate the origins of such non-linearities in transport, which is especially 
problematic since the system is necessarily driven far from equilibrium. We argue that for a wide 
class of systems the new ingredient that enters is the Schwinger mechanism — the production of 
carriers from the vacuum by the applied field — which is then balanced against a scattering rate 
which is itself set by the field. We show by explicit computation how this works for the case of the 
symmetric superfluid-Mott insulator transition of bosons. 

PACS numbers: 



The notion of quantum criticality provides one the few 
general approaches to the study of strongly correlated 
quantum many-body systems lj. The scale invariancc 
that characterizes the zero temperature critical point 
leads to characteristic universal power law dependences 
for various quantities in its proximity; these dependences 
can be computed within a continuum field theory 

While the power law dependences can generally be re- 
lated to an underlying set of (possibly unknown) expo- 
nents by scaling arguments, establishing the actual mech- 
anism that gives rise to them — in the sense of prescribing 
a controlled computation — is not always a trivial task. 
For example, much attention has focussed on understand- 
ing the real-time dynamics at non-zero temperatures 
where the textbook procedure of analytic continuation 
from Matsubara expressions typically yields incorrect an- 
swers and no insight. 

Here, we address another such question, that of trans- 
port at finite fields. Specifically, consider a quantum crit- 
ical point between an insulator and a metal or supercon- 
ductor/superfluid characterized by a correlation length 
that diverges as £ ~ 8~ v and an energy scale that van- 
ishes as A ~ 5 VZ , where S is a measure of the distance 
to the transition. Dimensional analysis implies that the 
zero-temperature conductivity obeys the scaling form 

a{S ,E) = E^(^ T ^j (1) 

in d spatial dimensions. Thence, generically, the sys- 
tem exhibits a non-linear current-voltage characteristic 
at criticality. The question of interest to us is how the 
system might produce such a response. Evidently, linear 
response theory or naive perturbation theory to higher 
orders is no help. Physically, the system must set up a 
steady state whose properties depend in singular fashion 
on E with no expectation that it resembles the steady 
state obtained in thermal equilibrium. 

We show that the properties of such weakly interact- 
ing fixed-points may be understood as follows: The ap- 
plication of an electric field leads to a biased growth 



of current carrying fluctuations by an analogue of the 
Schwinger mechanism by which an electric field pro- 
duces electron/hole pairs from the vacuum. This process 
is non-perturbative in the electric field. The fluctuations 
produced in this way scatter from one another due to 
the interactions at the fixed point, thus producing cur- 
rent relaxation which is again non-perturbative in the 
electric field at the fixed point. Together, these effects 
establish a steady-state distribution that carries a cur- 
rent. We implement this idea for the simplest quantum 
phase-transition with singular transport — the symmet- 
ric superfluid to Mott insulator transition of bosons in a 
periodic potential — and recover the dependence Q with 
z = 1. In particular, in d = 2 we find a linear conductiv- 
ity a = (7r/4)e* 2 //i, different from those calculated pre- 
viously |2j, |j] . The difference between these results arises 
due to the different regimes of frequency, temperature 
and electric field (uj,T,E) to which they apply and the dif- 
ferent physical processes important in each case. The lin- 
ear response calculated at zero temperature amounts 
to the limit (w -> 0, T = 0, E = 0), whilst the finite 
temperature, linear response|2( corresponds to (u> = 0, 
T =^ 0, E = 0). Our result follows from taking the fre- 
quency and temperature to zero at finite field; (u) = 0, 
T = 0, E 0). The non-commutativity of the limits 
lu — > and T — » was first recognized in0j. The re- 
sult presented here suggests that similar care should be 
taken with the electric field. Next, we outline our compu- 
tation which builds upon the important work by Damle 
and SachdevEJ on the finite-temperature transport of the 
same model 5J. 

Field Theory The critical region of the superfluid to 
Mott insulator transition with particle hole symmetry 
is described by a charged scalar field with a quartic 
interactionjlj, |2| : 

H = J d d x [tfn + V^V(/> + mW + A(0ty) 2 ] , (2) 

where <fi is the complex scalar field and Fl is its conju- 
gate momentum. These satisfy the usual commutation 
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relations [0(x, t), II(y, t)] = i<5(x — y). It is convenient to 
choose the bare interaction A to have its fixed point value 
u *j\3-d w jth momentum cutoff A [ij although we will not 
need the precise, regularization-dependent, value. At the 
zero-temperature critical point, the renormalized mass is 
zero which corresponds to a particular choice m* of the 
bare mass. The effects of applying an electric field, E, 
are included by minimally coupling to a vector potential; 

— * D(j> = (V + iA) 4>. We choose the gauge A = Et 
for initial convenience, later we switch to a scalar poten- 
tial. 

Mode Expansion The normal modes of this Hamil- 
tonian (in the absence of interaction) are charge den- 
sity fluctuations. These occur with positive and negative 
charges, corresponding to decrease or increase in charge 
density from the average. The first step in our analy- 
sis is to re-express (J2J in terms of creation and annihi- 
lation operators for these normal modes. The transfor- 
mation is standard for the Klein-Gordon theory. Using 
a) and a to represent the creation and annihilation of 
positively charged density fluctuations and b' and b for 
negatively charged fluctuations, the non-interacting part 
of the Hamiltonian reduces to 



dg(k,t) e k (t) 



Ho = 



A* (ot(k,t),6(-k,t)) 



( o(k,t) 
B k e k + B k J { &t(_k,t) 



(3) 



where e k = vF+n? is the mode energy and B k = 
((Ei) 2 + 2k.Et) /2e k . The Hamiltonian has an explicit 
time dependence arising from our choice of gauge used 
in coupling to the electric field. This time dependence 
is responsible for the production of fluctuations by the 
electric field. 

Schwinger Mechanism Let us ignore the interaction 
between the normal modes of the system. This will 
allow us to describe the pair production process, after 
which we will put back the interactions. Our first step 
is to diagonalize the Hamiltonian © using a Bogoliubov 
transformation. Because of the time dependence of the 
Hamiltonian, the transformation used to carry out this 
diagonalization must itself be time dependent. This has 
important consequences for the equations of motion. In 
the instantaneously diagonal basis, or adiabatic particle 
basis, the Heisenberg equations of motion for operators 
pick up extra terms from the time dependence of the 
Bogoliubov transformation matrix. We are concerned 
with the equations of motion for the regular and anoma- 
lous distribution functions; /(k,t) = (a'(k, i)a(k, t)) = 
(&t(-k,t)&(-k,t)) and g(k,t) = (b(-k, i)a(k, t)), where 
the angular brackets indicate averages over the state of 
system j^j. After transforming to the adiabatic particle 
basis, the equations of motion for these reduce to 



df(k,t) e k (t) 



dt 



2e k (i) 



(2/(k,t) + l)-2ie k (t)0(k,t) (4) 



where e k (t) = \J (e k + 5 k ) 2 - -B k = e k +Et is the mode 
energy in the adiabatic particle basis. These equations 
contain all of the ingredients necessary to describe the 
Schwinger mechanismQ. The terms proportional to 
e k (t) result from the time dependence of the Hamilto- 
nian and are responsible for pair creation. The second 
equation can be solved explicitly for g(k, t) and substi- 
tuted back into the first. The resulting equation contains 
a source term for the production of fluctuations, which 
includes oscillatory behavior coming from the quantum 
coherence of the pair production. If we ignore these 
transients, the source term may be replaced by a delta- 
function with the appropriate weight. The equation of 
motion for /(k, t) then reduces to 

d/(k,t) 



dt 



5(t + k.E/E 2 )e" 



-w(k± 2 +m 2 )/E 



(5) 



dt 



ek(t) 



Keg(k, t) 



The pair production described in this way may be un- 
derstood by analogy with Landau-Zener tunneling^] as 
follows: under the action of the electric field, the momen- 
tum of a charged excitation increases linearly in time as 
k + Et. The energy of the charged excitation changes in 
time accordingly; at large momenta it is proportional to 
|k + Ef|, and at a time t = — k.E/E 2 it passes though 
a minimum equal to m. This is very similar to the en- 
ergy dependence of modes in the Landau-Zener model. 
Except for the Bose enhancement factor 2/(k, t) + 1 the 
pair production equations (4) have precisely the same 
form as in Landau-Zener tunneling. 

Scattering and Steady State The steady creation of 
pairs from the vacuum will, in the absence of scattering, 
lead to a secular divergence of the current. Thus, consid- 
eration of the scattering is essential for understanding the 
steady state transport. Below d = 3, the critical behav- 
ior is controlled by the interacting, Wilson-Fisher fixed 
point. If we access the properties of this fixed point in 
a weak coupling expansion, such as the l/N expansion, 
which we use in this paper, then the leading order de- 
scription of the critical, steady state can be obtained by 
considering the scattering of the particles produced via 
the Schwinger mechanism. Heuristically, the pair pro- 
duction in (JjjJ leads to a growth in the current with a 
rate proportional to E~^~ , while scattering with a cou- 
pling of order 1 [N to N channels is expected (on dimen- 
sional grounds) to lead to a current relaxation rate of 
order (1/N)y/E; together these will reproduce with 
2 = 1. In the remainder of the paper we shall see how 
this insight can be turned into a computation within the 
quantum Boltzmann equation framework. 

Generally, 1/A^-expansions extend the number of 
modes of the model from its initial value to some large 
number N, allowing all of these modes to interact with 
one another. If N is taken to be very large, the interac- 
tion may be expanded perturbatively in l/N. Crucially, 



3 



in the present case, we couple the electric field to just two 
modes of our extended theory (i.e. one of the N/2 copies 
of the model) . The field induces fluctuations in these two 
modes via the Schwinger mechanism. These fluctuations 
scatter into the remaining N — 2 modes, thus allowing 
the modes that are coupled to the electric field to reach 
a steady state. We discuss moving beyond leading order 
later in the paper. The details of how this 1 /N expansion 
is set up are given in ||. 

The scattering integrals in our Boltzmann equations 
may be determined using several methods. In the case 
where only regular distribution functions are required, 
it is possible to use Fermi's Golden rule. Since we 
also have to consider anomalous distributions, this sim- 
plest approach does not work. It is, however, possible 
to determine the scattering integrals by using Heisen- 
berg's equations of motion with the 1/N Hamiltonian 
and time-dependent perturbation theory; this is after all 
how Fermi's Golden Rule is derived in the first place. 

1- Loop Scattering Integrals The 1-loop correction is 
particularly simple to calculate. After calculating the 
Heisenberg equations of motion including the interaction 
Hamiltonian and then using mean field theory, one finds 
that the only modification is to replace £?k with £>k + Sk, 
where the self-energy £k = m(E) 2 /(2ek)- An expres- 
sion for the 1-loop renormalized mass m(E) is given to 
lowest order in a 1/A^-expansion in Eq. JJJ below||. 
By making a new Bogoliubov transformation to the 1- 
loop renormalized adiabatic particle basis one may reduce 
the Heisenberg equations for the distribution functions 
to the same form as replacing the adiabatic mode 
energy ek(i) everywhere with the renormalized mode en- 
ergy, ek, m (t) = V ( k + Ei ) 2 + m(E) 2 . Notice that this 
modification points to a feedback between the production 
of fluctuations and the scattering between them. Since 
the rate of production depends exponentially upon the 
energy gap, it is suppressed by scattering — although its 
functional dependence on E is not changed. Of course, 
this scattering is not sufficient to establish a steady state 
and indeed, can be ignored at leading order in 1 /N alto- 
gether as will be clear below. 

2- Loop Scattering Integrals The next order scatter- 
ing integral is deduced from the Heisenberg equations of 
motion for the various mode operators followed by time 
dependent perturbation theory to lowest order in 1/N. 
To lowest order in 1/N and at zero-temperature, the 
dominant scattering is due to processes where the modes 
coupled to the E-field scatter off one another and into 
the N — 2 uncoupled modes of the extended theory. At 
this juncture, it is useful to make the change of variables 
t, q — > t, k, where k = q+Ei is the instantaneous momen- 
tum of the quasi-particle. Rather than label the normal 
modes by their initial momentum and keep the label fixed 
in time, we choose to label them with their instantaneous 
momentum. This is equivalent to changing gauge to a 
scalar potential. With this change of variables all explicit 



time dependence is removed, except for in the distribu- 
tion functions. In particular, e q (i) = £k = \/k 2 + m 2 . 
The resulting Boltzmann equations are 

[d t + -E.d k ]f(k,t) = +^^neg(k,t) 

+ 7k (t) [Jleg(k,t)- f(k,t)} 
[d t + E.d k ]g(k,t) = -2te k 2(k,i) + ^^^(2/(M) + l) 



£k 



- 7 k(i)/(k,i) 
with the supplementary definitions 



(6) 



7k (t) 



y/\k\+h f dk' 



N e k 



{2nf 



£k' 



-/(k',t) 



j . Wtt f dk f(k,t) 
m{E,t) = — / - ' ' +A 



N J (2tt) 2 e k 



(7) 



The damping factor 7(k, t) is derived allowing for the 
zero-temperature, critical propagation of the N — 2 un- 
coupled modes. 

Scaling and Current The above equations clearly per- 
mit a universal, steady state solution whose properties 
are governed only by the electric field and the fixed 
point value of the coupling. Writing the scaling forms 
f/g(k,E) = f/g(k/VE) = f/g(k), m = rhVE and 
7 = jVE) we see that the Eqs. © reduce to E inde- 
pendent equations for / /g(k). This by itself is sufficient 
to establish a current proportional to E d / 2 in dimensions 
d < 3. For d > 3 the fixed point is Gaussian and the 
necessity of including dangerously irrelevant scattering 
processes will modify this scaling. 

We can make progress on the actual solution by mak- 
ing two simplifications valid for the leading order (in 
1/N) computation of the current: we can ignore the 
mass renormalization and eliminate the second equation 
in favour of a local source term^(J- With these we find 
the greatly simplified and soluble equation 



E 



df(k) 
dk 



-7k/(k)+e-^/ B 5(k.E/£ 2 ) (8) 



with 7k and m(E) given by Eq. iJBJ in the case of steady- 
state distributions. 

The current carried in the steady-state is given by 



j = 2 



dk k\\ 



/(k) , 



(9) 



with an additional term involving g(k) being subdomi- 
nant in 1/N. Upon rescaling, this reduces to a form that, 
in 2-dimensions, is proportional to the electric field. The 
constant of proportionality is universal and may be calcu- 
lated in the 1/N limit to be a = (NTr/8)e* 2 /h 11], where 
N = 2 in the physical system. 
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Dissipation and higher orders The steady state that 
we have described involves a balance between pair pro- 
duction, their acceleration by the field and relaxation due 
to scattering. The latter processes need to relax the cur- 
rent, the number of charge carriers, as well as the energy. 
The first two require processes present in the Hamilto- 
nian but energy relaxation is present only at leading or- 
der where the infinitely many orthogonal modes act as 
a heat sink. In order that we be able to go to higher 
orders in l/N we need to understand how that problem 
is to be dealt with. Qualitatively we expect to mimic the 
logic of linear response theory, where Joule heating is an 
0(E 2 ) process which can be dealt with without invali- 
dating the O(E) result one computes. In our case the 
leading current response is 0(E d / 2 ) but the Joule heat- 
ing j • E is still down by a factor of E so the same logic 
is prima facie applicable. Explicitly, this can be done by 
constructing a solution that carries a heat current to the 
boundaries of the sample, e.g. transverse to the direction 
of electric current flow|l2j. Roughly, this can be thought 
of as a local "effective temperature" T c s(y) which drives 
a heat current nV 2 T e f{(y). By scaling, T c g cx \[E, Re- 
quiring that the variation in T e s be less than its mean 
value and factoring in the scaling form of k leads to the 
conclusion that such a solution can be constructed for a 
system with transverse dimension, L y , and electric fields 
satisfying L y ^/~E <C 1. We note that such a restriction is 
not unusual, a similar construction for a ordinary metal 
at finite temperature also yields a bound on system size 
and electric field^3|- While the estimate that we give 
would apply to a generic critical point with z = 1, our 
particular system is an even better bet for such a con- 
struction as k is infinite at all T due to an absence of 
energy current relaxation. 

Bose gases While this work was motivated by the 
fundamental question of understanding the non-linear 
quantum-critical states, we would be remiss if we did 
not note that the Mott transition described by (1) 
has been observed, remarkably enough, in atomic Bose- 
Einstein condensates placed in optical lattices of variable 
In this system, the role of the electric field 
is played by an intensity gradient in the optical field, or 
alternatively by an acceleration of the optical lattice. The 
universal current predicted within our analysis amounts 
to a steady flow of matter proportional to the acceleration 
of the optical lattice. Field gradients and accelerations 
in these systems are easily made large in the sense of our 
theoretical discussion. Whether the non-linear response 
discussed here can be observed given the complications 
of the confining potential and the differing requirements 
of equilibration in these systems appears to be a fit topic 
for further theoretical and experimental work. 
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